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Abstract. Bokermannohyla martinsi is an endemic anuran from the mountains of the southernmost portion of the Espi- 
nhaco range in southeastern Brazil. This region is known as the Iron Quadrangle (Quadrilatero Ferrifero) due to its iron- 
rich outcrops and is under intensive mining exploration, which threatens its caves. Here, we investigated whether males, 
females, and juveniles of B. martinsi use iron caves throughout the year, exhibit cave fidelity, distribute themselves differ- 
ently within caves, and prefer to use caves that are microclimatically stable. We sampled 10 caves during eight 15-day peri- 
ods regularly spaced throughout one year. Frogs used caves throughout the year and exhibited cave fidelity. They preferred 
the dysphotic zone, and young individuals remained closer to the entrance of the cave and to the ground, maybe due to 
their lower locomotory/climbing ability compared to adult frogs. Bokermannohyla martinsi preferred more stable caves 
(those with lower temperature and humidity variation), and the use of caves by females was more intensive during the dry 
and cold season, which may be related to their shorter permanence at breeding sites (streams). Our results show that caves 
are important habitats for B. martinsi, as it may be the case with other species, too, and their destruction should thus be 


evaluated and weighted with greater care when planning their mining. 
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Introduction 


A broad definition of “cave” as provided by CULVER & 
PIPAN (2019) is that it is “a natural opening in solid rock 
with areas of complete darkness, and larger than a few mil- 
limeters in diameter” Subsurface habitats include both 
large (actual caves) and small cavities (interstitial habitats; 
CULVER & PIPAN 2019). Cave use in animals can be classi- 
fied according to their dependence on a hypogean (subter- 
ranean) habitat. Based on an adaptation of the Schinner- 
Racovitza system (HOLSINGER & CULVER 1988) provided 
by SKET (2008), (1) troglobionts are obligatory inhabitants 
of the hypogean habitat; (2) eutroglophiles are epigean 
(surface) species with subterranean population(s); (3) sub- 
troglophiles may inhabit subterranean habitats for varied 
periods of time, but also need to exploit the epigean habitat 
for some resources; and (4) trogloxenes only occur spo- 
radically underground. However, varied combinations of 
morphological, physiological, and behavioural adaptations 
influence the level of association of species with caves, 


making it sometimes difficult to classify them in the ab- 
sence of complete information (HOWARTH & MOLDOVAN 
2018). 

Different from Europe, where the amphibian fauna in- 
cludes troglobiont species and other species well known 
to maintain populations in caves (e.g., BLACKBURN 2019, 
LunGut et al. 2020), many frogs in South America and oth- 
er regions have frequently been classified as accidental cave 
occupants when found in caves (TRAJANO 1986, TRAJANO & 
GNASPINI-NETTO 1990, PINTO-DA-ROCHA 1995, TRAJANO 
& BICHUETTE 2006). However, several studies have dem- 
onstrated that some species may use caves regularly dur- 
ing parts of or their whole lifecycle (RESETARITS JR. 1986; 
MILLER et al. 2008, ESPINO-DEL-CASTILLO et al. 2009, Su- 
WANNAPOOM et al., 2018), being thus classified as troglox- 
enes or subtroglophiles (sensu SKET 2008). 

Anurans were recorded at higher frequencies during au- 
tumn and winter in temperate caves, when environmental 
conditions are harsher outside, although they may also use 
cooler cave microhabitats during periods of high external 
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temperatures (RESETARITZ JR. 1986,FENOLIO et al. 2005, 
LUNGHI et al. 2017). However, standardized monitoring 
studies that relate anuran presence in caves to climatic var- 
iables are rare, especially from tropical regions (EsPINo- 
DEL-CASTILLO et al. 2009). Anurans are sensitive to tem- 
perature and humidity conditions due to their highly per- 
meable skin (WELLS 2007), and caves may be useful shel- 
ters for these animals during harsh weather (MATAVELLI et 
al. 2015). Caves have more stable microclimatic conditions, 
meaning they have little temperature variation and high 
humidity compared to the outside (CULVER & PIPAN 2019). 
There is great variation in microclimatic features among 
caves, with long caves with small entrances being more sta- 
ble than short caves with large entrances, as the latter are 
more influenced from outside conditions (FERREIRA 2004, 
PELLEGRINI et al. 2016). Although iron caves are usually 
not very deep (OLIVEIRA et al. 2011), they can have humid- 
ity and temperature conditions similar to those of caves 
from other lithologies, becoming relatively stable with in- 
creased distance from the entrance (FERREIRA 2005). 

Brazil hosts the highest amphibian richness in the world, 
with 1,091 known species (SEGALLA et al. 2019), 5% (54 spe- 
cies distributed in 19 genera and 11 families) of which have 
already been recorded from inside caves (MATAVELLI et al. 
2015). Iron caves have been used by 23 of these species (14 
genera and 10 families), more than any other cave lithol- 
ogy (MATAVELLI et al. 2015). Iron caves are usually shallow 
and relatively short, with a mean development of 25 m in 
the Iron Quadrangle in southeastern Brazil (OLIVEIRA et 
al. 2011). However, these caves shelter a high invertebrate 
diversity, including troglobiont species (FERREIRA 2005). 
The Iron Quadrangle is a unique tropical iron-rich rock 
outcrop with characteristic vegetation formations and high 
species diversity (JACOBI et al. 2007, LEITE et al. 2008), but 
is greatly endangered by the expansion of mining activities 
(Jacosi et al. 2007). 

Bokermannohyla martinsi (BOKERMANN, 1964) is en- 
demic to the Iron Quadrangle in southeastern Brazil (P1- 
NHEIRO et al. 2014), where it has been recorded only from 
above 800 m a.s.]. (CARAMASCHI et al. 2010). It breeds in 
permanent streams surrounded by riparian forests (Bo- 
KERMANN 1964, CARAMASCHI et al. 2010). Although it is 
listed as ‘of least concern’ in the IUCN Red List of Threat- 
ened Species, it is classified as ‘near threatened’ in the Bra- 
zilian List of Endangered Species (HADDAD et al. 2016) 
due to its restricted distribution and intensive habitat de- 
struction (mostly by mining activities) (CANELAS & BER- 
TOLUCI 2007, HADDAD et al. 2016). MATAVELLT et al. (2015) 
recorded B. martinsi in iron caves, and it has been record- 
ed also in limestone and quartzitic caves (M. C. M. An- 
drade, pers. obs.). The Iron Quadrangle is the region with 
the second highest number of iron caves in Brazil and the 
one with the highest concentration, encompassing 1083 
known caves in an area of 490 km’, following the Cara- 
jas region, with 1787 known caves in an area of 4196 km? 
(OLIVEIRA et al. 2011, RUBBIOLI et al. 2019). Information 
on the importance of these caves to B. martinsi, and other 
anuran species, could be useful to understand the general 


Table 1. Features of sampled caves at the Serra do Gandarela 
National Park, southeastern Brazil including cave identification 
(Cave id.), area, altitude above sea level (alt.), linear extension 
(LE), number of entrances, and summed extension of all en- 
trances (ent/extent), Environmental Stability Index (ESI), mean 
temperature amplitude (MTA), and mean humidity amplitude 
(MHA) measured within caves, and number of individuals of 
Bokermannohyla martinsi recorded from September 2017 to Au- 
gust 2018. 


Caveid. Area alt. LE ent/ ESI MTA MHA Number 
(m) (m) extent of frogs 
Apol-25 1 1607 12.3 3/3.26 1.41 0.87 1.21 3 
SG-16 1 1618 51.1 3/5.77 1.89 1.04 1.23 11 
SG-20 1 1601 33.3 5/6.2 1.04 1.23 2.83 3 
SG-51 1 1627 55.4 1/4.2 2.58 0.16 0.40 20 
SG-17 1 1602 7.4 3/4.75 0.29 1.67 4.22 0 
SG-18 1 1603 9.0 3/2.97 0.89 1.08 1.88 0 
SG-19 1 1602 8.2 4/5.27 0.90 0.75 1.48 
Gand-056 2 1517 9.9 1/2.75 1.77 0.61 3.32 5 
Gand-096 2 1554 45.0 5/5.18 1.91 0.85 2.06 20 
Gand-008 2 1552 98.9 8/6.38 2.66 0.50 0.94 39 


importance of caves as habitats for these animals, especial- 
ly in regions where both caves and species that use them 
are endangered. 

Here, we investigated whether the use of caves by 
B. martinsi is intentional and not occasional. Thus, we test- 
ed the predictions that individuals of B. martinsi (1) use 
iron caves throughout the year, (2) exhibit cave fidelity, 
(3) go deeper into caves and are less active during periods 
with harsher external climatic conditions, and (4) prefer 
to use caves that are more stable regarding temperature 
and humidity and depending on external climatic condi- 
tions. We accounted for potential differences between sex/ 
age groups (males, females and juveniles) in all analyses. 
We expected frogs to venture farther from the entrance in 
search of greater microclimatic stability when external hu- 
midity is low, which usually coincides with lower tempera- 
tures in the region. We also expected frogs to be less fre- 
quently active under such circumstances, as their breeding 
season (and expected higher activity levels) falls into the 
wet season. 


Materials and methods 
Study site 


We sampled 10 natural iron caves in two different are- 
as in the Serra do Gandarela National Park (20°673.24” S, 
43°39 48.33” W), southeastern Brazil (Fig. 1, Table 1). This 
region is covered by typical ironstone outcrop plant com- 
munities (JACOBI et al. 2007) with Atlantic Forest on the 
northeastern slope and Cerrado formations in the south- 
western portion (LAMOUNIER 2009). This region has a wet 
season from October through March and a dry season 
from April through September. 
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The sampled caves have no bodies of water inside them; 
their humidity is generated by rain that permeates through 
fissures in the rock. The two sampled areas are separated 
by about 1 km (Fig. 1). One lies higher up on the terrain 
and encompasses seven caves (1601-1627 m a.s.l.), with dis- 
tances between caves varying from 11 to 390 m. The other 
is on the slope and encompasses three caves (1517-1554 m 
a.s.l.), with distances between caves varying from nine to 
14 m. The linear extension of the 10 caves varies from 7.4 to 
98.9 m (mean = 33 m). 


Species sampling 


We sampled caves during eight periods between September 
2017 and August 2018, with a minimum interval of 30 days 
between consecutive sampling periods. Two 15-day sam- 
pling periods were undertaken within each of four three- 
month periods that we considered as early wet season (Oc- 
tober-December), late wet season (January-March), early 
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dry season (April-June), and late dry season (July-Septem- 
ber), respectively. During each 15-day sampling period, we 
sampled each cave for frogs on Days 1, 8 and 15, during 
the day. Sampling order was randomized among caves. We 
searched the whole extension of each cave from the floor to 
the ceiling visually with torchlights and listened for possi- 
ble frog calls. We only moved forward after checking every 
surface that we could see, and these encompassed almost 
the entire expanses of the caves. Sampling time per cave 
lasted from 3 min. (SG-17, SG-18, and SG-19) to 3 hours 
and 18 min. (GAND 008), with a mean of 24 min. Frog call 
recording and the manipulation of specimens are includ- 
ed in these times. The total sampling effort summed up to 
96 hours and 18 min/person. 

We measured the snout-vent-length (SVL) of each frog 
and identified the sexes of sexually mature individuals as 
follows: males are smaller than females and have hypertro- 
phied arms and an ossified prepolex (BOKERMANN 1964, 
Lima et al. 2013, MAGALHAES et al. 2018). Thus, we con- 
sidered as adult males those individuals showing these fea- 


638800 


® Sampled caves 
Zo National Park Serra do Gandarela 


— Brazil 
(J State limits 


Figure 1. Location of the two areas and the ten sampled caves in the National Park Serra do Gandarela, southeastern Brazil. The points 
in the maps with larger scales also represent the locations of the caves for reference (points overlap due to scale). 
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tures and as adult females those that were larger than the 
smallest adult male recorded (56 mm SVL) and did not 
present these features. Any individual smaller than 56 mm 
was considered juvenile. 

We photographed the dorsal region and limbs of eve- 
ry individual for subsequent identification, as their colour 
patterns are highly variable and facilitate individual recog- 
nition (Fig. 2). This photographic identification method 
has already been proved successful in other anuran species 
and renders unnecessary the use of invasive marking tech- 
niques (MIRANDA et al. 2005, KENYON et al. 2009, LAMA 
et al. 2011, Russo et al. 2018) and minimizes the handling 
of individuals (LUNGHI et al. 2020). We created a photo ar- 
chive for the identification of recaptured individuals with 
no need for additional manipulation. We released each 
frog at the exact same spot of capture after taking measure- 
ments and pictures. 

We recorded, for each frog, its distance from the nearest 
cave entrance (irrespective of its being usable to humans or 
not), height above the floor, and light incidence. We meas- 
ured the distances with a laser distance meter (Leica Dis- 
to™ D8) with a range from 0.05 to 200 m, and a precision 
of 1.5 mm. We classified light incidence into three catego- 
ries: (1) euphotic, close to the entrance with direct light inci- 
dence (that is, with insolation at some time during the day) 
and presence of photosynthetic organisms, (2) dysphotic, 
with indirect light incidence (with no direct insolation), and 


(3) aphotic, without any incidence of light (adapted from 
SoaRESs et al. 2013). As the caves are not long (Table 1), only 
the three longest actually had aphotic zones. We also record- 
ed the activity level of each individual at the moment of its 
discovery as inactive (resting with the eyes closed, body flat- 
tened to the surface, and limbs held close to the body, not 
responding to the torchlight) or active (awake with limbs 
extended to any extent and responsive to the torchlight). 


Environmental features inside and outside caves 


We recorded internal temperatures and humidity levels of 
the caves with data loggers (digital thermo hygrometers 
AKSO®, model AK174) placed on the floor in the middle 
of each cave’s linear extension, which corresponded to the 
area where most frogs were located, based on our previous 
observations. We also placed identical data loggers outside 
each of the two groups of caves to record temperature and 
humidity. We set them up on the first day of each 15-day 
sampling period and retrieved them on the last day. The 
data loggers recorded the respective values every 30 min- 
utes and we used mean values and standard deviations 
(SD) for each 15-day sampling period in our analyses (see 
Statistical analyses). We obtained data on mean tempera- 
ture, mean humidity, and mean daily temperature and hu- 
midity amplitudes. 


Figure 2. Individuals of Bokermannohyla martinsi illustrating the variation in dorsal colour patterns used for individual identification. 
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We calculated cave environmental stability using the 
Environmental Stability Index (ESI) proposed by FERREI- 
RA (2004). This index measures the influence of the exter- 
nal environment on the cave microclimate, taking into ac- 
count the dimensions of the cave, as well as the number, 
positions, and sizes of entrances (FERREIRA 2004, BENTO 
et al. 2016, PELLEGRINI et al. 2016). For caves with a single 
entrance, ESI = In (LE/EE), where LE = linear extension, 
and EE = entrance extension. The measurement of linear 
extension is obtained by summing up the longitudinal line- 
ar measurements of all conducts, irrespective of their incli- 
nation (CHABERT & WATSON 1981). We measured both LE 
and EE with the laser distance meter (Leica Disto™ D8). 

For caves with more than one entrance, we applied 


( (LE/ZEE) `\ 
ESI = In NEDER) 


where NE = number of entrances and DEE = mean of 
the distances from one entrance to all the others. Larger 
values of ESI indicate more stable caves. Caves with an ESI 
of between 1.00 and 3.50 are considered moderately stable, 
with large dysphotic zones (FERREIRA 2004). We compared 
cave ESIs with the mean values of temperature and humid- 
ity amplitudes (MTA and MHA, respectively, that corre- 
spond to mean standard deviations) measured by us inside 
the caves per 15-day sampling period in order to confirm 
the applicability of the index (following FERREIRA 2004). 
As we found ESIs to represent well cave stability concern- 
ing MTA (estimate = -0.77, t = -3.26, p = 0.011) and MHA 
(estimate = -0.44, t = -2.55, p = 0.034), we then used this 
index in subsequent analyses. 

We obtained rainfall data from the meteorologi- 
cal station at the Santuário do Caraça Private Reserve 
(20°5’53.6” S, 43°29'17.0” W). We used mean values of the 
ten days prior to each 15-day sampling period in our analy- 
ses (see Statistical analyses). We considered this period to 
adequately represent the rainfall before our sampling pe- 
riods, considering that it filled and maintained the water 
levels of the temporary streams where B. martinsi breeds at 
the study site (M. C. M. Andrade, pers. obs.) and thus likely 
influenced its activity levels. 


Statistical analyses 


We first confirmed individual frog identities through com- 
parison of the photographs taken and calculated the num- 
bers of recaptures, frog relocations between caves, and the 
numbers of individuals of each age/sex group in specific 
caves during each sampling period. 

We tested whether perch height above the ground and 
distance from the cave entrance (log-transformed) were 
influenced by external climatic factors. We used data from 
every record of cave use by individual frogs. We tested the 
roles of (1) sex/age group (adult male, female or juvenile), 
(2) activity level, (3) mean external temperature, (4) mean 
external humidity amplitude, and (5) mean rainfall in the 
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Table 2. Seasonal variation in cave use throughout the year rep- 
resented by numbers of males, females, juveniles, and all indi- 
viduals recorded in each seasonal period in all sampled caves at 
the National Park Serra do Gandarela, from September 2017 to 
August 2018. These numbers account for 429 records from 94 
individuals (several of them recorded in more than one seasonal 
period). 


ne ape dass Late dry Early wet Late wet Early dry 


(number of 

aera season season season season 
individuals) 

Total adults 44 28 24 26 
Adult males 26 19 15 15 
Adult females 18 9 9 11 
Juveniles 25 8 8 20 
Total individuals 69 36 32 46 


previous 10 days as explanatory variables for frog (1) dis- 
tance (log-transformed) from the nearest entrance to the 
cave, and (2) perch height above the ground as depend- 
ent variables. We did not include the mean external tem- 
perature amplitude, because it was correlated to the mean 
external humidity amplitude (Spearmans rho = 0.880, 
P < 0.001), nor mean external humidity, because it also 
correlated with the mean external humidity range (Spear- 
man’s rho = -0.905, p < 0.001). We used individual identity 
as a random variable to account for repeated encounters. 
We built Linear Mixed-Effects Models (LMM) using the 
package Ime4 (Bares et al. 2015) in R (R Core Team 2020). 
Model significance was assessed by comparison with a null 
model including only the random variable. 

We looked for cave features and external climatic fea- 
tures that influenced the use of caves by males, females and 
juveniles of B. martinsi throughout the year. We used data 
from each of the 10 caves in each one of the eight sampling 
periods, so that the dependent variable was the number of 
individuals (separate analyses for males, females, and juve- 
niles) detected in each cave during each sampling period. 
Repeated records of the same individual during the same 
sampling period were not considered. The explanatory var- 
iables were ESI, mean external temperature, mean exter- 
nal humidity amplitude, and mean rainfall during the pre- 
ceding 10 days. We used cave identity and sampling period 
as random variables. We built Generalized Linear Mixed 
Models (GLMM) with the function glmer (family Poisson, 
link = “log”) in the package Ime4 of R (BaTEs et al. 2015). 
We used the package MuMIn (BARTON 2019) to select the 
best models from the global model using AICc. We test- 
ed for residual overdispersion with the package DHARMa 
(HARTIG 2021). We used the package emmeans (LENTH 
2020) for post hoc tests. 


Results 


We obtained 429 records of 94 individuals of B. martinsi 
(Table 2) in seven out of the ten sampled caves. These in- 
cluded four caves from Area 1 and all the three caves of 
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Area 2 (Fig. 1), with a maximum of 39 individuals recorded 
in one cave (Table 1). The caves with more records and in- 
dividuals were the ones with greater linear extension (Ta- 
ble 1). Only seven out of the 94 recorded frogs used two 
different caves during the study. Two of these moved from 
cave SG-20 to cave SG-16 (Fig. 1), whose entrances are 
separated by 115 m (consecutive individual records were 
35 and 50 days apart). Three individuals moved between 
Gand-056 and Gand-o008, separated by 12 m (consecutive 
individual records were 5 to 58 days apart), and two oth- 
ers moved between Gand-o096 and Gand-008, separated by 
nine metres (both individual records were 31 days apart). 
We recorded 12 individuals (13%) during the four sam- 
pling periods, nine (10%) during three periods, 38 (40%) 
during two periods, and 35 (37%) during only one of the 
sampling periods. A mean of 59% of individuals were cap- 
tured during one period and again during the following 
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Figure 3. Number of individuals of Bokermannohyla martinsi 
captured during each season with the proportion representing 
recaptures from the previous season. 


Females 


period, with the greatest recapture rates coinciding with 
the beginning of the wet season (66%) and the beginning 
of the dry season (63%; Fig. 3). The numbers of new juve- 
niles increased at the end of the wet season and the begin- 
ning of the dry season, contributing to larger numbers of 
juveniles in caves recorded during these and subsequent 
months (February—August, Fig. 4), whereas the number of 
new adults increased from the beginning of the wet season 
until the beginning of the dry season, but we found no new 
adults at the end of the dry season. Individuals returned 
repeatedly to the very same spot where we had found them 
previously. 

The mean distance of individuals of B. martinsi from the 
cave entrance was 7.28 m, corresponding to the dysphotic 
zone. Compared to the null model, only sex/age group ex- 
plained variations in the distance of frogs to the cave en- 
trance (Chi-Square = 10.71, df = 2, p = 0.005) with juveniles 
remaining closer to the entrance than females (post-hoc 
estimate = 0.21, se = 0.06, t = 3.31, p = 0.004) and males, 
although this result only approaches significance (post- 
hoc estimate = 0.43, se = 0.06, t = 2.26, p = 0.067). Males 
and females did not differ in their distances to the cave en- 
trance (post-hoc estimate = 0.07, se = 0.05, t = 1.41, p = 
0.341; Fig. 5A). Again, only sex/age group explained vari- 
ations in the perch height of frogs above the ground (Chi- 
Square = 17.18, df = 2, p < 0.001) with juveniles remaining 
closer to the ground than both males (post-hoc estimate 
= -0.55, S€ = 0.15, t = 3.58, p = 0.002) and females (post- 
hoc estimate = -0.65, se = 0.16, t = 4.20, p < 0.001). Males 
and females did not differ in their perch height above the 
ground (post-hoc estimate = 0.10, se = 0.11, t = 0.93, p = 
0.625; Fig. 5B). No other models showed better explanatory 
power than the null model. 

The best model to explain cave use by males included 
only ESI (estimate = 2.95, se = 0.68, Z = 4.32, p < 0.001; 
Supplementary document 1). Models including rainfall 
and mean humidity amplitude had similar AICc values 
(Table 3), however neither rainfall nor mean humidity am- 
plitude alone explained a significant amount of variation 
in the number of males per cave. The best model to explain 


Juveniles 


Figure 4. Numbers of individual males, females, and juveniles of Bokermannohyla martinsi recorded per sampling period in ten sam- 
pled caves in the National Park Serra do Gandarela, southeastern Brazil. 
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the use of caves by females included ESI (estimate = 3.37, 
se = 0.96, Z = 3.51, p < 0.001) and mean external tempera- 
ture (estimate = -1.193, se = 0.07, Z = -2.81, p = 0.005) with 
females choosing more stable caves during colder periods 
(recorded during the dry season, Fig. 4). This model pre- 
sented deviant residual dispersion (Supplementary docu- 
ment 2). A model including just mean external tempera- 
ture (estimate = -0.193, se = 0.07, Z = -2.82, p = 0.004) had 
residuals fitted to expected values (Supplementary docu- 
ment 3). Models including rainfall and mean humidity 
range also had similar explanatory powers. Females tend- 
ed to use caves more often during drier periods, but mean 
humidity amplitude alone did not explain the cave use by 
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Figure 5. Distance (m; log-transformed) from cave entrance (A) 
and perch height (m) above the ground (B) of males, females 
and juveniles of Bokermannohyla martinsi recorded in ten sam- 
pled caves in the National Park Serra do Gandarela, southeastern 
Brazil. 
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females. The best model to explain the use of caves by juve- 
niles included just ESI (estimate = 1.36, se = 0.51, Z = 2.70, 
p = 0.007; Supplementary document 4). Other models in- 
cluding rainfall, mean temperature, and mean humidity 
amplitude had similar explanatory powers, however nei- 
ther of these variables alone explained a significant amount 
of the variation in the cave use by juveniles. 


Discussion 


We recorded an intensive and continuous use of caves by 
the endemic frog Bokermannohyla martinsi, with seven out 
of the ten sampled caves hosting frogs (four of those dur- 
ing all sampling periods throughout the year) and 63% of 
all individuals recorded being found in caves during more 
than one sampling period. Site fidelity was also observed, 
with 52 (88%) recaptured frogs being found in the same 
caves. Recruitment of juveniles (70%) occurred at the end 
of the dry season, and they were likely from the previous 
rainy season, as some of them still had tail remnants. 

Even though we do not know which proportion of the 
whole population uses caves, at least 94 individuals oc- 
curred in this habitat, which is more than we usually find 
in breeding habitats (streams) during year-round sam- 
plings in other localities (AFONSO & ETEROVICK 2007, P. C. 
Eterovick, pers. obs.). Thus, the high abundance of B. mar- 
tinsi in the studied caves corroborates the hypothesis that 
they use these habitats as shelters. During the whole study 
we never saw any potential predator inside the caves. 

Bokermannohyla martinsi used more intensively the 
dysphotic zone, which has high humidity and moder- 
ate temperatures (Biswas 2014) as is usually observed in 
troglophilic and trogloxenic species (LUNGHI et al. 2015, 
2017), including other anurans (OSEEN & WASSERSUG 
2002, LUNGHI et al. 2018). In these sectors, the frogs may 
also benefit from decreased visibility to visually orientated 
predators (DUELLMAN & TRUEB 1994, OSEEN & WASSER- 
SUG 2002) and also from moister and more stable microcli- 
matic conditions compared to the euphotic zone (HETEM 
et al. 2012, EVERALL et al. 2015, LuNGuiet al. 2018). Adult 
individuals of B. martinsi were found to have moved far- 
ther into the caves than juveniles, which could be due to 
the former's better locomotory performance resulting from 
their larger sizes and greater familiarity with the cave habi- 
tat due to previous visits, enabling them to reach more sta- 
ble sectors. 

Within the caves, juveniles tended to be closer to the 
ground than adults of B. martinsi. Spots higher up may be 
chosen during periods of lower activity, as has been ob- 
served in Eleutherodactylus longipes (BAIRD, 1859) (EsPI- 
NO-DEL-CASTILLO et al. 2009), however, active and inac- 
tive individuals of B. martinsi did not differ regarding their 
perch heights above the ground. For their part, juveniles 
may position themselves closer to the ground due to their 
less developed climbing skills, as opposed to adult frogs. 
Morphological features of B. martinsi such as digital disks, 
a camouflage colour pattern, and nocturnal habits have also 


Caves as frog habitats 


Table 3. Models built to evaluate the influence of cave ESI (Environmental Stability Index) and external climatic features: rainfall, 
external mean temperature (extmeantemp), and external mean humidity amplitude (extMHA) on use of caves (number of individuals 
using caves per sampling period) by males, females and juveniles of Bokermannohyla martinsi at the National Park Serra do Gandarela, 
southeastern Brazil, from August 2017 to August 2018. The selected model is presented in boldface. All models that differed from the 
best model by less than two AICc units are presented. Significance values for the increased explanatory power of the selected model 
against the corresponding null model are presented (see text for analysis details). 


Models to explain number of males in caves per sampling period df  loglik AICc Delta Weight p 
males ~ ESI+(1|cave)+(1|sampling) 4 -81.72 171.97 0.00 0.38  <0.001 
males ~ ESI+rainfall+(1|cave)+(1|sampling) 5 -81.22 173.25 1.29 0.20 

males ~ ESI+extMHA+(1|cave)+(1|sampling) 5 -81.41 173.64 1.67 0.16 

Models to explain number of females in caves per sampling period 

females ~ ESI+extmeantemp+(1|cave)+(1|sampling) 5 -67.59 146.00 0.00 0.30  <0.001 
females ~ ESI+rainfall+extMHA+(1|cave)+(1|sampling) 6 -66.76 146.68 0.68 0.21 

females ~ ESI+rainfall+(1|cave)+(1|sampling) 5 -68.07 146.94 0.95 0.19 

females ~ ESI+rainfall+extmeantemp+(1|cave)+(1|sampling) 6 -67.03 147.22 1.22 0.16 

females ~ ESI+rainfall+extmeantemp+extMHA+(1|cave)+(1|sampling) 7 -66.00 147.55 1.55 0.14 

Models to explain number of juveniles in caves per sampling period 

juveniles ~ ESI+(1|cave)+(1|sampling) 4 -64.75 138.04 0.00 0.24 0.007 
juveniles ~ ESI+extmeantemp+extMHA(1|cave)+(1|sampling) 6 -62.77 138.70 0.66 0.17 

juveniles ~ ESI+extMHA(1|cave)+(1|sampling) 5 -64.21 139.22 1.19 0.13 

juveniles ~ ESI+extmeantemp (1|cave)+(1|sampling) 5 -64.25 139.32 1.28 0.13 

juveniles ~ ESI+rainfall+extMHA(1|cave)+(1|sampling) 6 -63.11 139.38 1.34 0.12 


been found in other anuran species that use caves (BISWAS 
2014, LURiA-MANZANO & RAMiREZ-BAUTISTA 2017) and 
may represent important adaptations for cave users. 

Although the positioning of frogs within the caves did 
not depend on external climatic conditions, the number of 
frogs using caves was higher during the dry season, when 
a greater recruitment of juveniles is expected and exter- 
nal climatic conditions are harsher. During this period, 
both adult frogs and juveniles may seek shelter in caves. 
This effect was significant for females of B. martinsi, which 
used caves more often during colder and dry months. Al- 
though the model including ESI to explain the use of caves 
by females had overdispersed residuals, this is not unusu- 
al in GLMMs and does not mean the model is unusable 
(HARTIG 2021). This seasonal pattern was also observed 
for Eleutherodactylus longipes in a cave in Mexico (Esp1- 
NO-DEL-CASTILLO et al. 2009). In the Crimean mountains 
(Ukraine), amphibian species also used caves non-acciden- 
tally during the winter (non-breeding season), although 
breeding-related activities are likely to occur within caves 
in some species (TURBANOV et al. 2019). 

The lower numbers of records of B. martinsi in caves 
during the rainy summer is likely due to the frogs being 
outside in search of breeding habitats, with males remain- 
ing at calling sites for longer periods than females (PCE 
pers. obs.). Frogs may then seek temporary shelter in spots 
closer to streams during the daylight hours between con- 
secutive nights of breeding activity. Although B. mar- 


tinsi uses permanent streams for reproduction (BOKER- 
MANN 1964, CARAMASCHI et al. 2010), we found tadpoles 
in ephemeral streams in the vicinities of the studied caves. 
In Area 1, the closest streams are 120-180 m from the en- 
trances of the caves, whereas in Area 2, the closest stream is 
about 160 m distant. These distances between streams and 
caves are within the possible reach of these amphibians, as 
is suggested by their moving 159-290 m between breed- 
ing habitats and terrestrial habitats (SEMLITSCH & BODIE 
2003). However, we recorded B. martinsi also in caves with 
perennial bodies of water (one quartzitic and three dolo- 
mitic caves in the Serra do Gandarela and one quartzitic 
cave in the neighbouring Serra do Caraga). Thus, the pos- 
sible use of caves for reproduction by B. martinsi is worth 
investigating. 

Bokermannohyla martinsi used caves more intensive- 
ly during dry periods. Although B. martinsi can breed all 
year round (with peaks in the rainy season) in permanent 
streams (ETEROVICK & BARATA 2006), its breeding activi- 
ties at our study site are restricted to the rainy season, when 
temporary streams carry water for tadpole development. 

Bokermannohyla martinsi did not use caves with the 
smallest (< 1.00) ESI values and used more frequently 
caves with greater environmental stability. The preference 
of amphibians for caves with greater humidity and envi- 
ronmental stability was also observed in other field stud- 
ies (LUNGHI et al. 2014) and experimentally (LUDDECKE 
2003). Amphibian susceptibility to water loss via the skin 
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(WELLS 2007) makes the use of caves advantageous due to 
their high humidity and stability (Biswas 2014), facilitat- 
ing water condensation and absorption through the skin, 
even at sites with low water availability (Tracy et al. 2011). 
Additionally, prey availability and lower predation risk are 
likely advantages of caves (BISWAS 2014). 

This study underscores the importance of caves as habi- 
tats for anurans, corroborating many others that have dem- 
onstrated these animals not to be just occasional cave vis- 
itors (RESETARITS JR. 1986, FENOLIO et al. 2005, ESPINO- 
DEL-CASTILLO et al. 2009, BISWAS 2014, LUNGHI et al. 2014, 
MATAVELLI et al. 2015, TURBANOVet al. 2019). The practical 
application of such knowledge in the process of environ- 
mental licensing should afford caves greater appreciation 
and acknowledge their need of preservation, which is es- 
pecially necessary in the Iron Quadrangle, considered one 
of the major areas of iron mining in the world (PINHEIRO 
et al. 2014), where many caves exist and are threatened with 
destruction through mining. Brazilian legislation (Execu- 
tive order #6640, 7 November 2008) already attributes a 
higher conservation value to caves hosting troglobionts 
and/or that represent essential habitats for populations of 
threatened species, amongst other features. The current as- 
sessment of the importance of cave habitats to the species 
that use them is, however, not satisfactory, making it dif- 
ficult to apply the law properly (Cruz & P1L6 2019). We 
hope our results, by adding to the knowledge of the impor- 
tance of caves as habitats for species not usually considered 
cave inhabitants, will contribute to a greater appreciation 
of caves in the conservation of biodiversity. 
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Supplementary data 


The following data are available online: 


Supplementary document 1. Residual diagnostics for the model 
selected to explain the use of caves by males of Bokermannohyla 
martinsi in ten sampled caves in the National Park Serra do Gan- 
darela, southeastern Brazil. 
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Supplementary document 3. Residual diagnostics for a model in- 
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the National Park Serra do Gandarela, southeastern Brazil. 


Supplementary document 4. Residual diagnostics for the model 
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